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In borocarbide Erffel^C the phase transition to the commensurate spin density wave at 2.3 K leaves 1/20 
part of Ising-like Er spins practically free. Vortices polarize these spins nonuniformly and repolarize them when 
moving. At a low spin relaxation rate and at low bias currents vortices carrying magnetic polarization clouds 
become polaron-like and their velocities are determined by the effective drag coefficient which is significantly 
bigger than the Bardeen-Stephen (BS) one. As current increases, at a critical current J c vortices release polar- 
ization clouds and the velocity as well as the voltage in the I-V characteristics jump to values corresponding to 
the BS drag coefficient. The nonuniform components of the magnetic field and magnetization drop as veloc- 
ity increases resulting in weaker polarization and discontinuous dynamic dissociation depinning transition. As 
current decreases, on the way back, vortices are retrapped by polarization clouds at the current J r < J c , 

PACS numbers: 74.25.Wx, 74.70.Dd, 74.25.Ha 



The family of quaternary nickel borocarbides (RE)Ni2B2C 
(RE is rare earth magnetic ion) is an interesting class of crys- 
tals which exhibits both superconductivity and magnetic order 
at low temperatures. flJO. A number of the crystals in that 
family develop antiferromagnetic order below the Neel tem- 
perature TV below the superconducting critical temperature 
T c . It has been recognized time ago that superconductivity 
coexists quite peacefully with the antiferromagnetic order as 
the spatial periodicity of magnetic moments is well below the 
superconducting correlation length. In contrast, the ferromag- 
netic order, antagonistic to the Cooper pairing, leads to dra- 
matic changes in both magnetic and superconducting orders 
in the coexistence phase, for a review see Ref. Q]. That is 
why interest in the compound ErNi2B2C with T c = 1 1 K and 
Tfi — 6 K peaked when it was realized that below the phase 
transition from incommensurate spin density wave (SDW) to 
commensurate SDW at T* = 2.3 K the phase with weak fer- 
romagnetic ordering may emerge. [5,61 From neutron scat- 
tering data it was concluded that in ErNiaBaC below the 
incommensurate SDW develops with effective Ising spins ori- 
ented along the a-axis and with the wave vector Q = 0.5526/?*, 
where b* - 2n/b and b is the lattice period along the b-axis. 
Q [H At T* the transition to the commensurate phase with 
Q = 0.55b* leaves one out of 20 spins free of SDW molecular 
field. These Er spins with the magnetic moment fi = 7.8/kb are 
easily polarizable by the magnetic field along the a direction. 
The spin magnetization in the magnetic field H = 2000 G in 
temperature interval 2 K - 4 K follows M sp /H sx /uM s /(k B T), 
where M s as 56 G, see Fig. 4 in Ref. [6 |. This value, M s = ym, 
corresponds to magnetization when all "free" spins with the 
concentration n order ferromagnetic ally, The same value M s 
was obtained by extrapolation of the magnetization at temper- 
ature 2 K in fields H > 1500 G to H -> 0. @ Note, that 
the Hall probe measurements below T* without an applied 
field found magnetic internal field much lower than M s and 
no spontaneous vortex lattice was seen. [ 10 1 High polarizabil- 
ity of spin system in ErNi2B2C is a key point for our following 
discussion. 



As hope to observe remarkable consequences of weak fer- 
romagnetic phase coexisting with superconductivity waned, a 
few puzzles on ErNi2B2C behavior at low temperatures re- 
mained. First, it was discovered by measuring the hystere- 
sis in the M - H loops and transport measurements that new 
pinning mechanism develops below T* for which the critical 
current increases as temperature lowers. ]9][TT] Second, neu- 
tron scattering data in applied magnetic field H close to H || c 
have shown that vortices deviate randomly from the direction 
of the magnetic field inside the crystal. [12| The vortex de- 
viations increase proportionally to l/T as temperature drops 
from 4 K to 1 .6 K and is nearly independent of the magnetic 
field in contrast to the usual behavior when effect of disorder 
drops with field. So far, no explanation of surprising temper- 
ature and field dependence of the critical current and disorder 
has been offered. 

To explain these data we propose a new mechanism of pin- 
ning - formation of polaron-like vortices dressed by the po- 
larization cloud of magnetic moments. The polaronic mecha- 
nism is inherent to all magnetic superconductors but it is best 
pronounced when the magnetic system is highly polarizable, 
as in the case of ErNiaBaC below 4 K. To clarify this mech- 
anism, we recall that the magnetic field is nonuniform within 
the vortex lattice being strongest near the vortex cores. Con- 
sequently, the polarization of the magnetic moments is also 
nonuniform. When vortices move they should repolarize the 
magnetic system, otherwise they would lose the energy gained 
by polarization (the Zeeman energy). The process of repolar- 
ization depends on the dynamics of magnetic system. In the 
following we consider the relaxation dynamics of free spins 
in ErNi2B2C. The repolarization process is controlled by the 
relaxation time t which should be compared with the charac- 
teristic time a/v needed to shift the vortex lattice moving with 
the velocity v by the vortex lattice period a = (<Do/B)'^ 2 (<t>o 
the flux quantum and B the magnetic induction). For r » a/v 
the magnetic moments slow down strongly the vortex motion. 
At some critical velocity and critical current, the vortices are 
stripped off polarization clouds. The corresponding jump in 
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FIG. 1 . (color online) Schematic view of the vortex lattice in the pres- 
ence of free Ising magnetic moments along the a axis. The vortex lat- 
tice is tilted from the applied magnetic fields in the ac plane due to the 
polarization of the magnetic moments. The vertical columns show 
the vortex cores. The polarized magnetic moments are nonuniform 
in space due to the spatial modulation of the vortex lattice magnetic 
field. Due to the Lorentz force F L vortices move along the x axis. In 
moving lattice, there is a phase shift between the magnetic induction 
(dashed line) associated with the vortex lattice and the magnetization 
(solid line) caused by the retardation in the response of magnetic mo- 
ments to the vortex magnetic field. 



velocity is strongly pronounced for a large t's. As current 
decreases, the vortices become retrapped again at the current 
J r < J c . Since the voltage V cc v, the I-V characteristics show 
hysteresis. The physics here is similar to that of a polaron 
with vortices playing the role of electrons and the magnetic 
polarization in place of phonons lfl3l . 

The E1-N12B2C crystals have orthorhombic structure below 
7V with domains where a- and b-axes change by 90° in neigh- 
bouring domains. We consider first clean single-domain crys- 
tal and later will discuss the effect of domain walls. We con- 
sider the vortex lattice induced by applied magnetic field H 
tilted by the angle 8 with respect to the crystal c axis. We 
choose the z axis along the direction of vortex lines at rest and 
x axis in the aoplane, see Fig. [T] The vortex line deviates 
from the applied field H due to the magnetic moments fl2l . 
In static situation the direction of vortex lines is determined 
by the effective field H + 4nM. Here M is the spatial aver- 
age of the magnetization. We denote by a the angle between 
vortex lines and the c axis. The deviation of vortex lines from 
the applied magnetic field was also discussed in Ref. Ifl4ll in 
the case of spontaneous ferromagnetic order of the magnetic 
moments. 

In the London approximation the magnetic field of the vor- 
tex lattice inside the crystal is (r = x, y) 



B 



cos(G • r) 
A 2 G 2 + 1 ' 



(1) 



where G are reciprocal vectors of the square lattice, A is 
the superconducting penetration length renormalized by the 
magnetic moments and B is the averaged magnetic induction. 



Here we ignore anisotropy of the penetration length. As re- 
vealed by neutron scattering, vortices form square lattice in 
ErNi 2 B 2 C. d 

In the Lagrangian the interaction between vortex lines at 
R, = (Xi,yi), and the magnetic moments is determined by the 
term 



£mt{R/,M) 



drB : (Ri - r, t)M z (r, t), (2) 



where we describe the magnetic moments in the continuous 
approximation via the magnetization M z (r, t), because dis- 
tance between free spins, 35 A, 1 8 1 is much smaller than the 
London penetration length A, about 500 A. Ifl2ll We ignore the 
pair breaking effect [15| of the magnetic moments because 
they suppress Cooper pairing uniformly as distance between 
free spins is much smaller than the coherence length, and thus 
moments do not introduce pinning. We also neglect the ef- 
fect of disorder in crystal lattice. The equation of motion for 
vortex lines is 



dRi <9£ w {Ri,R;} d£ int {Ri,M} 



dt 



dR, 



OR: 
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where r\ is the Bardeen-Stephen drag coefficient, £ W {R,, R y } 
is the vortex-vortex interaction, next term describes the force 
acting on the vortex line R, from the magnetic moments and 
Fl = Q>oJ/c is the Lorentz force due to the bias current J. 
The force due to magnetic moments is the same for all lines 
and the vortex lattice moves as a whole. The motion of vortex 
lattice center of mass, u(t), along the jt-axis is described by 
the equation 



du d 
dt du 



J dvB z (x + u, 



y, f)M z (r, t) 



+ Ft 



(4) 



Using linear response approach to relate magnetization with 
the magnetic field we obtain 



du d 



>1 



dt du 
B z (r',t') + F L 



J drdr'B z (x + u, y, t) J dt'x zz (r - r', t - f) x 



(5) 



Here ^ zz (r, f) is the dynamic susceptibility of the magnetic 
moments. The vortex lattice moves with a constant velocity, 
// = vt, in the steady state t » r. Integrating over coordinates 
and time we obtain 



rjv • 
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where ;f z; (k, u) is the dynamic magnetic susceptibility in the 
Fourier representation. We see that the magnetic moments 
affect strongly vortex motion if a) the resonance condition 
v • G = Q(G) is fulfilled, where £2(k) is the frequency of mag- 
netic excitations with the momentum k and Q(k) » F(k), 
where T(k) is the relaxation rate of excitation, and b) dynam- 
ics of magnetic system is dominated by relaxation, Q(k) < 
F(k), favouring formation of the polaron. In the former case, 
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discussed in Ref. [ 16 1, the magnetic moments renormalize the 
vortex viscosity at high velocities when alternating magnetic 
field of vortices is able to excite magnons. Here we consider 
the latter case of free moments described by the relaxation 
dynamics with^-^k, to) given by 



1 



1 — iu)T 



X = 



(7) 



at temperatures T below 4 K for ErNi2B2C (effect of oredered 
spins will be discussed below). 

We introduce dimensionless quantities by expressing t in 
units of t and u in units Gl 1 = a/(2n). In the summation over 
G y in Eq. (|6ji, we account only for the dominant term with 
G y = 0. For GoA » 1 we find equation for the dimensionless 
vortex velocity 



fjv = F(v) + F L , F(v) = -v[n 2 /3 + v 2 - 7TVCOth(7r/v)], 
where we introduced the dimensionless parameters 
4ttV 4 * Ajt 2 GqA a J 



il = 



Ft = 



^cOosin a 



(8) 



(9) 



The asymptotic behaviour of F(v) is F(v) as -v at v « 1 
and F(v) « -n 4 /(45v) at v » 1. The electric field due to 
the motion of the vortex lattice is given by E — Bv/c, and 
we obtain the I-V curves from Eq. (|8j, as depicted in Fig. [2] 
As v increases, the v — J curve changes from a weak current 
dependence v = F L /(1 + fj) to a stronger and usual Bardeen- 
Stephen behaviour v = Fi/fj. In Fig. [2]at fj = 0.1, the I-V 
curve is hysteretic. Upon ramping up the bias current, the sys- 
tem jumps to the usual Bardeen-Stephen (BS) Ohmic curve at 
a current J c , where electric field increases discontinuously by 
the factor 1 / 77 at 77 <k 1. The jump, identifying experimen- 
tally as depinning transition, is caused by the dissociation of 
the vortex-magnon polaron. It is very similar to the dissoci- 
ation of usual electron-phonon polaron in high electric fields 
as described theoretically [ 17] and confirmed experimentally 
in metal oxides 1 18 ]. Upon decreasing the current the vortices 
are retrapped by the polarization clouds at a threshold current 
J r and the vortex lattice moves with a significantly enhanced 
viscosity at lower currents. 

The critical current J c and retrapping current J r follow from 
the equation for the velocity rj-dF{v)/dv = 0. The maximum 
of dF(v)/dv is 0.297, thus hysteresis exists for fj < 0.297. The 
calculated J c , J r and corresponding electric fields are shown 
in Fig. [3] At small fj the critical current is 



0.03 



Xc<S>o sin a 
GnA A 



(10) 



J c decreases with temperature as J c ~ l/T and decreases with 
the magnetic field as J c ~ 1 / yB. 

We note, that above T*, in the incommensurate SDW, some 
spins experience quite weak SDW molecular field. Thus, they 
are polarized by vortices and exhibit polaronic effect and pin- 
ning. This explains increase of pinning in ErNi2B2C as T 
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FIG. 2. Calculated I-V curves for fj = 0.1 and r) - 0.5. For f] = 
0.1, the system shows hysteresis in the I-V curve while for fj = 0.5, 
no hysteresis is present. The green dotted line denotes the unstable 
solution. 



decreases below 7V, see Ref. 0, as well as pinning in the 
holmium borocarbide below T N . lfP9l 

Let us consider the origin of the jumps at J c and J r . The 
dependence of the magnetization on the velocity of moving 
vortices is 



M z (r, v, t) = xB sin 2 a V 



cos[G-r-jg(v)] 
4^ W 2 G 2 + l)[l+(G. v vr) 2 ] 3 / 2 



(11) 



with tan(/3) — G x vt. Nonuniform component of the magne- 
tization and thus the polarization effect decrease with veloc- 
ity. On the other hand, the retardation between the magnetic 
field and the magnetization, as described by the phase shift 
B(v), increases with the velocity. This positive feedback and 
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FIG. 3. Dependence of the critical current J c and retrapping current 
J r , and corresponding electric fields E c and E r on fj. 
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the increase of retardation with velocity ensure discontinuous 
transitions at J c and J r . 

Strong pinning due to the polaron mechanism requires 
small parameter fj. It is expressed via r as 77 as 
10~ u s/(Tsin 2 a), where we have used the BS drag coefficient 
iIbs — ^n/i^^^Pn) with the coherence length £ « 13 nm 
|[T2l and the normal resistivity p„ = 5 /ifi-cm at T c . EUll The 
relaxation time t in ErN^I^C is long because the dynamics 
of majority of spins is strongly suppressed by the formation 
of the SDW molecular field as was found by the Mossbauer 
measurements. EH The relaxation time drops very fast below 
10 K and reaches the value t » 5 x 1(T 10 s at T = 5 K, but 
data at lower temperatures were not reported. Thus the only 
information we have so far is 77 < 0.02/ sin 2 a. 

The critical current for ErNi2B2C reported in Ref. [9| for 
B = 0.1 T, T = 2 K is about 250 A/cm 2 . Eq. ([ToJ» gives 
such current value at a — 2.5°. In experiment the applied 
magnetic field was close to the c axis, but the precise angle 9 
was not reported. [9 |. The estimate of the order 1° is reason- 
able, but quantitative comparison is not convincing as we do 
not know r and thus 77 below 2.3 K. Hence, the real check of 
polaronic mechanism should be by measuring the I-V charac- 
teristics. We predict hysteretic behavior in ErNi2B2C, strong 
dependence of voltage and of the critical current on the angle 
9, at least for 9 » 0.15. Note, that the critical current reaches 
values as high as 10 6 A/cm 2 at high angles at T — IK and 
B = 0.1 T. 

The effect of ordered spins on the vortex motion is similar 
to that described in Ref. 16 for an antiferromagnet. When 
Cherenkov condition v ■ G > O(G) is met, excitation of 
magnons results in enhanced viscosity 77. This occurs at high 
velocities, due to a gap in the magnon spectrum, and thus at 
high currents J > J c , leading to a voltage drop in comparison 
with the BS result. 

Let us discuss now the effect of disorder on the vortex lines 
direction observed in magnetic fields tilted with respect to the 
c-axis. [ 12 1 Due to domain structure the Ising spins are polar- 
ized only in half domains where vortex lines follow direction 
of the effective field H + 4;tM, while in others they are along 
H. The random change of angle of vortex directions with re- 
spect to average angle is 2nM/B = 2nx sin a. It increases as 
1 IT, when T drops, in agreement with the results of Ref. lfT2l 
and data for M sv /H mentioned above. When vortices cross a 
domain wall between domains with different magnetization, 
they need to repolarize magnetic moments at currents below 
the critical one. This slows additionally vortex motion, but 
such an effect is smaller than that accounted for previously be- 
cause domain size is much bigger than the distance between 
the magnetic moments. Domain walls also cause pinning of 
vortices as seen in Bitter decoration patterns. [22 1 Importantly, 
for a nonzero 9 small part of vortices experience domain wall 
pinning and this part drops with 9. In contrast, polarization 
effect increases with 9 and this helps to separate domain wall 
pinning from the polarization one. 

In conclusions, vortices in magnetic superconductors polar- 
ize magnetic moments and become dressed and polaron-like. 



At low currents and long spin relaxation time the nonuniform 
polarization induced by vortices slows their motion at currents 
for which pinning by crystal lattice disorder becomes ineffec- 
tive. As current increases above the critical one, vortices re- 
lease nonuniform part of the polarization and the velocity as 
well as the voltage in the I-V characteristics jump to much 
higher values. At decreasing current vortices are retrapped by 
polarized magnetic moments at the retrapping current which 
is smaller than the critical one. The results of such polaronic 
mechanism are in qualitative agreement with the experimental 
data (9j [T2| but measurements of the I-V characteristics are 
needed to establish the quantitative agreement and confirm 
validity of such a model for Er borocarbide. The polaronic 
mechanism should be at play also in Gd and Tb borocarbides 
superconductors in the commensurate SDW phase. It may be 
present in Tm borocarbide above T^ and in cuprate supercon- 
ductors (RE)Ba2Cu3C>7, where magnetic RE ions positioned 
between superconducting layers interact weakly with super- 
conducting electrons and order at very low Neel temperatures 
of the order 1 K. HI 
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